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ABSTRACT 

We present an analysis of the growth of black holes through accretion and bulges through star formation in 
33 galaxies at the centers of cooling flows. Most of these systems show evidence of cavities in the intracluster 
medium (ICM) inflated by radio jets emanating from their active galactic nuclei (AGN). We present a new 
and extensive analysis of X-ray cavities in these systems. We find that AGN are energetically able to balance 
radiative losses (cooling) from the ICM in more than half of our sample. Using a subsample of 17 systems, 
we examine the relationship between cooling and star formation. We find that the star formation rates are 
approaching or are comparable to X-ray and far UV limits on the rates of gas condensation onto the central 
galaxy. The remaining radiative losses could be offset by AGN feedback. The vast gulf between radiative losses 
and the sink of cooling material, which has been the primary objection to cooling flows, has narrowed and, in 
some cases, is no longer a serious issue. Using the cavity (jet) powers, we place strong lower limits on the 
rate of growth of supermassive black holes in central galaxies, and we find that they are growing at an average 
rate of ~ 0.1 M Q yr -1 , with some systems growing as quickly as ~ 1 M yr" 1 . We find a trend between bulge 
growth (star formation) and black hole growth that is approximately in accordance with the slope of the local 
(Magorrian) relation between black hole and bulge mass. However, the large scatter in the trend suggests that 
bulges and black holes do not always grow in lock step. With the exception of the rapidly accreting supercavity 
systems (e.g, MS 0735.6+7421), the black holes are accreting well below their Eddington rates. Most systems 
could be powered by Bondi accretion from the hot ICM, provided the central gas density increases into the 
Bondi radius as p oc r~ l . However, if the slope of the gas density profile flattens into a core, as observed in 
M87, Bondi accretion is unlikely to be driving the most powerful outbursts. 

Subject headings: cooling flows — galaxies: active — galaxies: clusters: general — X-rays: galaxies — 
X-rays: galaxies: clusters 



1. INTRODUCTION 

The intracluster gas at the center of a majority of 
galaxy clusters has a cooli ng time les s than 10 10 yr 
llEdge. Stewart. & Fabianlll992tlPeres etafl[i 998). In the ab- 
sence of a source of heat, this gas should cool, resulting in 
a slow inwar d flow of material know as a "cooling flow" 
(Fabianl [l994l) . While early observations of this gas made in 
X-rays supported this picture, observations of the clusters at 
other wavelengths did not. Optical data implied star forma- 
tion rates in the central g alaxy of only a few percent of th e de- 
rived cooling rates (e.g., McNamara & O'Connell 1989), and 
radio obser vations found less cold gas than predicted (e.g., 
lEdgdl20"0lh . However, recent high-resolution X-ray spectra 
from XMM-Newton do not show the features e xpected if large 
amounts of gas are coolin g below kT ~ 2 ke V ( Peters on et alJ 
120031 RCaastra et alJ2004l) . In addition, high spatial resolution 
Chandra X-ray Observatory images reveal that AGN can have 
a large heating effect (via jets and winds) on the ICM, sup- 
plying enough heat in some systems to offset radiation losses 
( Birza net alJl200 4). The emerging picture of cooling flows 

1 On leave from the School of Engineering Physics, University of Wollon- 
gong, Wollongong, NSW 2522, Australia 



is one in which most (but not all) of the cooling is roughly 
balanced by heating from AGN feedback, resultin g in a mod- 
erate cooling flow (e.g.fPedlar etalJl!^90llSokereta!]l2001t 
lBinne\l200l ISoke J20o3T ~ 

In this regulated-cooling scenario, net cooling from 
the ICM would lead to condensation of gas onto the 
central galaxy, driving the star formation observed in 
many systems (e.g.. Ijonhstone. Fabian. & Nulsenl 119871 
McNamara & 0'Connellll989l) . If this scenario is correct, the 
star formation rates should on average be comparable to the 
rate of gas observed to be condensing out of the ICM. Studies 
of a small number of systems with reliable star formation and 
cooling rates have shown that the rates are converging, and in 
some cases are in rough agreement (e.g., |M cNamaraj|200: 
iMcNamara. Wise. & Murravl l2004t iMcNamara et alJ 1200^ 
The quality and quantity of data from Chandra and XMM- 
Newton now makes it possible to construct a significantly 
larger sample of such systems to better understand the pos- 
sible connection between star formation and net cooling in 
cooling flow clusters. 

Additionally, the high-resolution data from Chandra 
are useful in studies of the nature of the feedback mecha- 
nism that may be preventing large amounts of intracluster 
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gas from cooling. Chandra images of galaxy clusters 
have revealed many large-scale interactions between 
the ICM and the central AGN, the bes t-known exam- 
ples of which are t he Per se us cluster ( Bohringer_etaL| 
1993} iFabian et all l2000t ISchmid t. Fabian . & Sander J 
2002; Fabian. Wilman. & Crawford 2002; Fabian et al. 
2003a b, 200 61), Abell 20 5 2 felanton et alJ 1200 U 
Blanton. Sarazin. & McNamaral 120031), an d Hydra A 
(McNamar aetalJ 120001 David et alT 1200 It iNulsen et all 
120021 12005bl) . In these systems, the radio jets of the AGN 
have pushed out cavities in the cluster's atmosphere, creating 
surface-brightness depressions in X-ray images that are 
correlated with the lobes' radio emission, such that the radio 
emission fills the depression in X-rays. The lower emis- 
sivities of the depressions imply that they are low-density 
cavities in the ICM, a nd therefore should rise buoyantly in the 
cluster's atmosphere JChurazov et al J 120011) . By measuring 
the surrounding pressure and volume of the cavities using 
the X-ray data, one can derive the work done by the radio 
source on the ICM in inflating the cavities, giving a direct 
measurement of the non-radiative energy released during 
the outburst. Measurements of this energy, combined with 
measurements of the star formation and cooling rates, can 
be used to investigate possible feedback scenarios that may 
govern the growth of the central dominant galaxy (CDG, as 
distinct from the more strictly defined cD) and its central 
supermassive black hole. 

Such feedback has implications for the more general prob- 
lem of galaxy formation. The large-scale distribution of mass 
in the Universe is well modeled by the standard hierarchical 
cold dark matter (CDM) cosmology JWhite & Reesll978l) . In 
this model, larger dark matter halos form through the merg- 
ing of smaller ones, while their gravitationally bound baryons 
cool and c ondense into the progenitors of the galaxies we 
see today JCoIe1fl99l lBlanchard. Valls-Gabaud. & Mamon 
1992). This picture successfully explains much of the ob- 
served matter distribution. However, a persistent problem in 
simulations that include only gravitational heating is a failure 
to reproduce the truncation of the h igh-luminosity end of the 
galaxy luminosity function (Benso n et al.ll20 03). This prob- 
lem stems from excessive cooling of baryons in the cores 
of halos, resulting in a population of massive galaxies, far 
larg er than even the enormous cD galaxies at the cores of clus- 
ters (Siiacki & Springe] 2006). Instead of residing in the cD 
galaxy as predicted by simulations, most baryons are found in 
the hot ICM. 

This problem may have a solution in non-gravitational 
heating by supernovae and AGN. Supern ovae are es sential 
for enriching the ICM to observed levels (Metzler & Evrard 
1994; Borgani etl^PzOOl and may play a significant heat- 
ing role in smaller galaxies, but their feedback ener- 
gies are too smal l and localized to truncate cooling in 
massive galaxies (Borgani et al. 2002). Furthermore, in 
the closely related preheating problem, they have diffi- 
culty supplying enough heat to boo st the entropy level of 
the hot gas to the observed levels (Wu, Fabian, & Nulsen 
1200(1 iBorgani et all 120051) . Energetically, AGN heating ap- 
pears to be the most likely mechanism to severely reduce 
the supply of gas from the hot ICM in galaxies above 
a certain size and to explain observed entropy p rofiles 
(Bens on et alJl2003n ScannaDieco & Oh 2004; Donahu e et al.l 
20061 1 Voit & Donahuel2005HVoit. Kav. & Brvanl2005l) . 

Current theory posits that AGN are powered by the accre- 
tion of material onto a central black hole. Gravitational bind- 



ing energy of the accreting material powers radiation and out- 
flows from AGN as the black hole grows. The relativistic jets 
that are revealed by their synchrotron emission are a prod- 
uct of this process. The remaining accreting material goes 
to increasing the mass of the black hole. In a sense, AGN 
are the "smoking guns" of black hole growth. The fraction 
of accreted power that re-emerges from an AGN and its parti- 
tioning between radiation and outflows is not well unde rstood , 
but probably depends on accretion rate (e.g.,|R ees et a|J1222l 
Naravan & Yill994t lAbramo wicz et alll995tlChurazov et alJ 
2005). We can place lower limits on the AGN's power using 
estimates of the power required to create the cavities associ- 
ated with the radio lobes. This power may then be used to 
infer the minimum growt h rate of the black ho le. 

As present e d by lFerrarese~ & Merritt (2000) and 
Gebhardt et al. (200(J), a correlation exists between the 
mass of the central black hole (Mbh) and the velocity disper- 
sion (cr) of the galaxy's bulge. This correlation suggests that 
the large-scale properties of the galaxy and the small-scale 
prope rties of the black hole are related (the "Magorrian rela- 
tion", [Magorrian et a ll!998l) . Estimates of the current growth 
rates of the black hole may be compared to the large-scale 
properties of the galaxy (such as the star formation rate) to 
trace the present day impact of bulge and black hole growth 
on this connection. 

In this paper, we use star formation rates, ICM cooling 
rates, and AGN heating rates for a sample of cooling flows to 
investigate the relationships between star formation and cool- 
ing and between the growth rates of black holes and their host 
galaxies. We assume Hq = 70 km s -1 Mpc" 1 , Qa = 0.7, and 
Om = 0.3 throughout. 

2. THE SAMPLE 

Our sample was drawn primarily from the iBirzan et alJ 
(2004) sample of cooling flows whose central galaxies show 
evidence of AGN activity as revealed by cavities in X-ray im- 
ages. We have supplemented this sample with a number of re- 
cently discovered cavity systems and one non-cavity system 
(A1068) with high-quality star formation data. In total, our 
sample comprises 3 1 CDGs, 1 group dominant galaxy (HCG 
62), and 1 giant elliptical (M84). Table^hsts the properties 
of the sample and references for publications that discuss the 
cavities or X-ray data. All the systems in our sample were ob- 
served with Chandra and have data publicly available in the 
Chandra Data Archive. The sample ranges in redshift from 
Z = 0.0035 to z = 0.545 and varies in its composition from 
groups to rich clusters. We note that our sample is biased in 
favor of cavity systems; therefore, conclusions drawn from 
this sample may not apply to cooling flows as a whole. 

3. DATA ANALYSIS 

The following section describes the reduction and analysis 
of data used in this paper. Briefly, the cavities were identified 
and their sizes and projected radial distances were measured 
from Chandra X-ray data. The temperature and density of 
the ICM as a function of radius, used to find the pressures at 
the projected locations of the cavities, were also derived from 
Chandra data. The internal pressure of the cavities was de- 
rived under the assumption that the cavities are approximately 
in pressure balance with the surrounding medium. The cav- 
ity's pressure, size, and position were then used to find the 
mean cavity power. Black hole growth rates were inferred 
from the cavity power under the assumption that accretion 
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TABLE 1 
Sample Properties. 



<T c a M bu i ge c 



System 


7 


CDG 


(km s -1 ) 


M K b 


M R b 


(10 11 M e ) 


References 


A85 


0.055 


PGC 002501 


340±9 


-26.72 ± 0.04 


-24.80 ±0.08 


31±1 


10, 26 


A133 


0.060 


ESO 541-013 




-26.36 ±0.06 


-24. 18 ±0.05 


17.9±0.4 


18, 19 


A262 


0.016 


NGC 708 


255±8 


-25.65 ±0.03 


-22.77 ±0.02 


4.9±0.1 


5 


Perseus 


0.018 


NGC 1275 


247±10 


-26.23 ± 0.04 


-24.25 ±0.01 


19.2±0.1 


12, 13, 15, 43 


2A 0335+096 


0.035 


PGC 013424 




-26.15 ±0.05 


-24.18 ±0.13 


18±1 


31 


A478 


0.081 


PGC 014685 




-26.64 ±0.07 


-24.66 ±0.10 


28±1 


44, 48 


MS 0735.6+7421 


0.216 


PGC 2760958 




-26.37 ±0.17 


-24.51 ±0.10 


24±1 


36 


PKS 0745-191 


0.103 


PGC 021813 




-26.82 ±0.09 


-24.63 ±0.10 


27±1 


21 


4C 55.16 


0.242 


PGC 2506893 




-26.10±0.13 


-24.75 ±0.50 


30±8 


22 


Hydra A 


0.055 


PGC 026269 


322±20 


-25.91 ±0.06 


-24.67 ±0.05 


28.2±0.7 


9, 33, 38, 40 


RBS 797 


0.350 












45 


Zw 2701 


0.214 


PGC 2401970 




-26.26 ±0.17 


-24.75 ±0.10 


30±1 


2 


Zw 3146 


0.291 


2MASX J 10233960+041 1116 




-27.67 ±0.14 






2, 23 


A1068 d 


0.138 


PGC 093944 




-26.71 ±0.08 


-25.07 ±0.21 


41±4 


35,50 


M84 


0.0035 


M84 


298±2 


-24.69 ±0.02 


-22.62 


4.3 


16 


M87 


0.0042 


M87 


341 ±3 


-25.55 ±0.02 


-23.61 


11 


16,51 


Centaurus 


0.011 


NGC 4696 


257±6 


-26.02 ±0.02 


-23.70 ±0.01 


11.6=1=0.1 


42, 45 


HCG 62 


0.014 


NGC 4778 




-25.26 ±0.03 






49 


A1795 


0.063 


PGC 049005 


294±10 


-26.50 ±0.08 


-23.86 ±0.10 


13.4±0.6 


11 


A1835 


0.253 


2MASX J14010204+0252423 




-27.36 ±0.14 






37, 46 


PKS 1404-267 


0.022 


IC 4374 


258±7 


-25.30 ±0.03 


-22.93 ±0.20 


5.7±0.5 


25 


MACS J1423.8+2404 


0.545 












1 


A2029 


0.077 


PGC 054167 


366±9 


-27.44 ±0.05 


-24.39 ±0.02 


21.9±0.2 


7 


A2052 


0.035 


UGC 09799 


259±11 


-26.27 ±0.06 


-23.62 


11 


3,4 


MKW 3S 


0.045 


NGC 5920 




-25.55 ±0.06 


-23.67 ±0.05 


11.2=1=0.3 


30, 32 


A2199 


0.030 


NGC 6166 


302±4 


-26.37 ±0.03 


-24.03 ±0.03 


15.7±0.2 


24 


Hercules A 


0.154 


PGC 059117 




-26.45 ±0.11 


-23.95 ±0.50 


15=1=4 


39 


3C388 


0.092 


PGC 062332 


365±23 


-26.24 ±0.06 


-24.46 ±0.50 


23±6 


28, 29 


3C401 


0.201 


PGC 2605547 






-23.43 ±0.50 


9±2 


41 


Cygnus A 


0.056 


PGC 063932 




-26.70 ±0.06 


-23.47 ±0.35 


9±2 


27, 47 


Sersic 159/03 


0.058 


ESO 291-009 




-26.26 ±0.10 


-23.68 ±0.39 


11±2 


52 


A2597 


0.085 


PGC 071390 


222±18 


-25.55±0.11 


-23.49 ±0.21 


9±1 


8, 34 


A4059 


0.048 


ESO 349-010 


296±49 


-26.74 ± 0.05 


-25.00 ±0.02 


38.2±0.4 


5, 19 



REFERENCES. — (1) Allen et al. 2004; (2) Bauer et al. 2005; (3) Blanton et al. 2001; (4) Blanton et al. 2003; (5) Blanton et al. 2004; (6) Choi et al. 2004; (7) 
Clarke et al. 2004; (8) Clarke et al. 2005; (9) David et al. 2001; (10) Durret et al. 2005; (11) Ettori et al. 2002; (12) Fabian et al. 2000; (13) Fabian et al. 2003a; (14) 
Fabian et al. 2005; (15) Fabian et al. 2006; (16) Finoguenov & Jones 2001; (17) Forman et al. 2005; (18) Fujita et al. 2002; (19) Fujita et al. 2004; (20) Heinz et al. 
2002; (21) Hicks et al. 2002; (22) Iwasawa et al. 2001; (23) Jeltema et al. 2005; (24) Johnstone et al. 2002; (25) Johnstone et al. 2005; (26) Kempner et al. 2002; (27) 
Kino & Kawakatu 2005; (28) Kraft et al. 2006; (29) Leahy & Grizani 2001; (30) Mazzotta et al. 2002; (31) Mazzotta et al. 2003; (32) Mazzotta et al. 2004; (33) 
McNamara et al. 2000; (34) McNamara et al. 2001; (35) McNamara et al. 2004; (36) McNamara et al. 2005; (37) McNamara et al. 2006; (38) Nulsen et al. 2002; 
(39) Nulsen et al. 2005a; (40) Nulsen et al. 2005b; (41) Reynolds et al. 2005; (42) Sanders & Fabian 2002; (43) Sanders et al. 2005; (44) Sanderson et al. 2005; (45) 
Schindler et al. 2001; (46) Schmidt et al. 2001; (47) Smith et al. 2002; (48) Sun et al. 2003; (49) Vrtilek et al. 2002; (50) Wise et al. 2004; (51) Young et al. 2002; 
(52) Zakamska & Narayan 2003. 

a Central stellar velocity dispersions were taken from the HyperLeda database; when more than one measurement was available, a weighted average was used. For 
the purposes of the buoyancy-age calculation, when no velocity dispersion was available, the average value for our sample ((er) = 295 km s _1 ) was adopted. 

^Total magnitudes from the 2MASS catalog (S'-band) or HyperLeda catalog (W-band), corrected for Galactic extinction, /^-correction, and evolution (see text for 
details). 

c Bulge mass calculated from the R-band absolute magnitude. Errors reflect uncertainties in Mr only. 

^The Chandra image of A1068 does not show evidence of cavities. A1068 is included because of the large starburst in the central galaxy. 

onto the central black hole fuels the outburst. Lastly, high- 
quality star formation and cooling rates were taken from the 
literature. We also use lower-quality cooling rates derived 
from Chandra data. Unless otherwise noted, errors and up- 
per limits are la values. 

3.1. Chandra X-ray analysis 

All systems were observed with the Chandra ACIS detector 
and data were obtained from the Chandra Data Archive. The 
data were reprocessed with CIAO 3.3 using CALDB 3.2.0 and 
were corrected for known time-dependent gain and charge 
transfer inefficiency problems. Blank-sky background files, 
normalized to the count rate of the source image in the 10-12 
keV band, were used for background subtraction. 2 

3.1.1. Temperature and Density Profiles 
2 See http://asc.harvard.edu/contrib/maxim/acisbg/ 



Spectra were extracted in elliptical annuli centered on the 
X-ray centroid of the cluster with eccentricity and position an- 
gle set to the average values of the cluster isophotes. Weighted 
response files were made using the CIAO tools MKWARF 
and MKACISRMF or MKRMF (MKACISRMF was used for 
all observations taken at the -120 C focal plane temperature; 
MKRMF was used for all other observations). 

Gas temperatures and densities were found by depro- 
jecting the spectra with a single-temperature plasma model 
(MEKAL) with a foreground absorption model (WABS) us- 
ing the PROJCT mixing model in XSPEC 11.3.2, between 
energies of 0.5 keV and 7.0 keV. The redshift was fixed to the 
value given in Tabled an d the foreground hydro gen column 
density was fixed to the Galactic value o flDickev & Lockmanl 
( 1990), except in the case of 2A 0335+096 and A478, when 
a significantly different value was required by the fit. In these 
two cases, the column density in each annulus was allowed to 
vary. 
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TABLE 2 

Cavity and ICM Properties Derived from Chandra Data. 



^cav,tot Lx M coo i c Z-cool '"cool 

System (10 58 erg) (10 42 ergs-') (10 42 ergs-') (M Q yr 1 ) (lO^ergs" 1 ) (kpc) 



A85 

A133 

A262 



1 2+'- 2 
24+J 1 

o 1 3+010 



37 +37 

620+1° 
9.7+Z-5 



365 ± 20 
106 ±2 
1 1 l+°- 4 



I8+1 3 
5±3 
<0.7 



Perseus 


19+ 20 

-5 


150+10° 

-30 


554 ±2 


20+| 


2A 0335+096 


l.l+io 

-0.3 


24+f 


338 ±2 


29+ 7 


A478 


1.5+^ 


100+ 8 o 

-20 


1440 ±10 


40+ 40 , 

—20 


MS 0735.6+7421 


1600+™° 

-600 


6900+'™ 

-2600 


450 ±10 


20+ 2 o* 

-10 


PKS 0745-191 


-10 


1700+i 4 o° 

-300 


2300 ±30 


170 ±90 


4C 55.16 


12+] 2 


420+^0 

-160 


640 ±20 


70 ±30 


Hydra A 


64+?? 
-11 


430+ 2 o° 

-50 


282 ±2 


16±5 


RBS 797 


38+ 50 


1 200+' 700 


■> 1 nn+'oo 

Jluu -130 


?no+ 490 * 
zuu_ 180 


Zw 2701 


350+«o 

-200 


6000+ 8 ?oo 

-3500 


430+?o 

-30 


< 8* 


Zw 3146 


-1 10 


5800+6800 

-1500 


3010^n 


590+!»o 

-170 


A 1068 




20 e 




< 48 


M84 


0.003+0°^ 

-0.002 


-0.6 


0.07 ±0.01 


0.038 ±0.00 


M87 


o n?n+° 014 

U - UZU -0.003 


6 0+ 4 2 


8 10+ 003 
5 -' u -0.04 


1 2+ 01 


Centaurus 


0.060+0 05' 

""-0.015 


7.4+5-8 

—1.8 


28.1 ±0.3 


2.7+0-2 

-O.l 


HCG 62 


0.046+00" 

-0.028 


3.9+f > 


1.8±0.2 


< 0.3 


A 1795 


4 7+6.6 
+•'-1.6 


160+23° 


625+6, 


8+1 3 


A1835 


47+ 50 
+ '-16 


1800+$° 


3160+1° 




PKS 1404-267 


12+0 .15 


20+J 6 


27 ±1 


5±2 


MACS J1423.8+2404 


9Q+52 


1400+2500 
1 -90O 


2290 ±30 


140+no 

1HU -90 


A2029 


4 8+ 2 - 7 
+'°-0.1 


S7+ 49 


1160±10 


< 1.9 


A2052 


1 7+ 2 - 3 


150+2°° 


97 ±1 


7 0+09 


MKW3S 


38+f 


410+ 420 


104 ±2 


5+ 3 

J -2 


A2199 


7 c+6.6 
'■ J -1.5 


2 7 o+250 


142 ±1 


<3 


Hercules A 


31+ 40 


310^0 


210+' 


< 58* 


3C 388 


5 2+ 7 - 5 

-2.1 


200+ 2 80 


27+J 


< 3* 


3C401 


11+20 


650+' 200 

OJ -420 


37+f 


12+5* 


Cygnus A 


84+™ 


1300!™ 


420 ±4 


31« 


Sersic 159/03 


25+f 


780+go 


220 ±6 


15±9 


A2597 


3.6+^f 


u '-29 


470+8 ? 




A4059 


3 0+ 25 


96+89 

'-35 


93±1 


3+2 



30+2° 
3±2 
< 0.3 
21+8 
13±4 
405° 
12+' 3 
230 ±120 

70 ±20 
13±4 

250+ 400 
<6 

680+;™ 



01 2+000 3 

u - ulz -o.ooi 
1 l+oi 

4.3 ±0.2 

<0.1 

10i 14 



3±1 
90+70 

< 3 

■x 4+0.5 
•'■+-0.2 

5+ 3 

J -2 

<3 
<46 

<2 
7±3 
50 ±10 
9±5 
30+ 30 
2±1 



142 

93 

57 

90* 
135 
150 
141 
176 
162 
109 
185 
135 
186 
152 

10 

26* 

54* 

33 
135 
156 

83 
187 
140 

87 
120 

91 
104 

55 

62 

91 
136 
128 

85 



a Cavity power calculated assuming 4pV of energy per cavity and the buoyancy timescale. 
I'Bolometric luminosity between 0.001 and 100 keV inside r c00 i. 

c Net cooling rate to low temperatures. Values marked with an asterisk where derived from observations with a low number of counts 
inside the cooling radius (< 15000) and are therefore less reliable. 

^Radius of the cooling region, inside which the cooling time is less than 7.7 X 10 9 yr (except for values marked with an asterisk, 
which correspond to the radius at the chip's edge). 

c For A1068, the cavity power was calculated from the u = 1400 MHz radio power as Pca^tot ~ 1500 X ^mHz^V-- 



3. 1 .2. Cavity Power 

Cavities seen in the X-ray emission of clusters allow a di- 
rect measurement of the non-radiative energy output via jets 
from the AGN. This measurement is independent of the radio 
properties and is the most reliable available, since its deriva- 
tion rests on only a few, well-understood quantities. The 
radio-emitting jets are understood to be displacing the ICM at 
the location of the cavities, doing work against the surround- 
ing plasma, as well as supplying thermal energy to the radio 
plasma that fills the lobes. The total energy required to create 
a cavity is equal to its enthalpy, given by 

E cav = —^-pV, (1) 
(7-1) 

where p is the pressure of the gas surrounding the cavity, V is 
the cavity's volume, and 7 is the ratio of specific heats of the 



gas inside the cavity. For a relativistic gas, 7 = 4/3, and the 
enthalpy is 4pV. We assume this value of 7 for all subsequent 
calculations involving E cav . 

The cavity's size and position were m easured following the 
procedures used in Birzan et al. (2004). The cavity was as- 
sumed to be in pressure equilibrium with its surroundings, 
and hence its pressure was taken to be the azimuthally av- 
eraged value at the projected radius of its center. The cav- 
ity's age was estimated in three ways: by assuming the cavity 
to be a buoyant bubble that rises at its terminal velocity, by 
assuming that the bubble moves outward at the local sound 
speed, and by assuming the cavity's age is governed by the 
time required for material to refill the volume of the cavity as 
it moves outward (for a detailed description of our analysis, 
see Birzan et al. 2004). These ages typically differ by factors 
of 2-4. For simplicity, we use the buoyancy age as the esti- 
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mate of cavity's age. This estimate is probably an upper limit 
on the age of the cavity (neglecting projection effects), since 
the cavity is expected to move outward supersonically during 
the early, momentum-dominated phase of the jet. 

The mean jet power required to create a cavity or cavity pair 
is then 

Pcav=^f, (2) 

where t is the average time between outbursts. This time 
is known only for a few objects, such as Perseus, for 
which the interval between outbursts may be estimated from 
the presence of mu ltiple generations of cavities and ripples 
(Fabia net all200 6). In objects with only a single set of cavi- 
ties, which make up most of our sample, the cavity's buoyancy 
age is used for t. 

As noted above, the buoyancy age is likely to be an overes- 
timate of the true age, and time evolution in the output of the 
AGN's jets can lead to underestimates of the amount of to- 
tal energy traced by the cavities. Furthermore, the discovery 
of shocks in a number of deep X-ray images of clusters (e.g., 
Cygnus A, Wilson et al. 2003; NGC 4636, Jones et al. 2002; 
MS0735+7421, McNamara et al. 2005; Hercules A, Nulsen 
et al. 2005a), which typically represent a comparable amount 
of energy as that contained in the cavities, may mean that cav- 
ities trace a fraction (~ 50%) of the energy of a typical out- 
burst. Lastly, our cavity powers do not include the radiative 
luminosity of the AGN. Therefore, our estimates of P cav rep- 
resent a lower limit to the total power of the AGN. 

Table [2] lists the total cavity energies and the asso ciated 
powers for the systems in our sample (see Table IA1I in the 
Appendix for the properties of each cavity). For A1068, 
in which no cavities are apparent in the X-ray image, the 
outburst power was estimated using the v = 1400 MHz ra- 
dio flux from the NVSS survey (Shoo = 23.1 ± 1.1 mJy), as 
^cav.tot ~ 1500 x i^hA, the average relation found from the 
sample of Bir zan et all J2004I) for radio-filled cavities. 

3.1.3. X-ray and Cooling Luminosities 

As in iBirzan et alJ J2004I) . we wish to compare the cavity 
powers to the heating rates required to balance losses from the 
ICM due to X-ray emission. These losses may be estimated 
as the difference between the total X-ray luminosity and the 
luminosity of gas cooling to low temperatures (i.e. out of the 
X-ray band). In this analysis, we define the cooling radius as 
the radius (or semi-major axis if elliptical annuli were used) 
within which the gas has a cooling time less than 7.7 x 10 9 
yr (the time since z = 1 , representative of the time that the 
cluster has been relaxed and a cooling flow could become es- 
tablished). For those systems in which the cooling radius lies 
beyond the chip's edge, we use the radius at the chip's edge 
as the cooling radius. 

The deprojection described in Section l3.1.1l was performed 
again, and the bolometric flux of the MEKAL component in- 
side the cooling radius was used to calculate the X-ray lumi- 
nosity, Lx- The same model, with the addition of a cooling- 
flow component (MKCFLOW), was used to obtain an esti- 
mate of the net cooling rate and the associated cooling lumi- 
nosity (Lcooi) of gas cooling to low temperatures (found by fix- 
ing the MKCFLOW low temperature to 0.1 keV). In the case 
of A1835, the spectra were of insufficient quality to obtain a 
reliable cooling rate (see McNamara et al. 2006). Table|2]lists 
the luminosities and cooling rates derived from Chandra data. 
We use these rates, and those from XMM-Newton and FUSE 



(described in Section l3~7l . as estimates of the net cooling rate 
of the ICM, which we com pare to the star formation rate of 
the central galaxy in Section l4~3l 

3.2. Black Hole Growth Rates 

The energies and ages described in Section 13.1.21 may be 
used to infer the minimum growth rate of the black hole as- 
suming the cavities were created by AGN jets fueled by ac- 
cretion onto the central black hole. Although the luminous 
energy radiated by the AGN is not included in the cavity en- 
ergies and must also be fuel ed by accretion, the se systems are 
radiatively inefficient (e.g., Birzan et al. 2004), and the con- 
tribution of radiation to the current total power is negligible. 
The outbursts might pass through a radiatively efficient phase 
during their initial stages, but this phase could not have been 
long-lived, since cluster AGN do not now show the quasar- 
like activity, and should not therefore affect our results signif- 
icantly. We stress that our black hole growth rates are lower 
limits; any energy in excess of the jet energy would result in 
underestimates of the average accretion rates, but we expect 
this effect to be small. 

The jets are produced through the partial conversion (with 
efficiency e) of the gravitational binding energy of the accret- 
ing material into outburst energy. The energy required to cre- 
ate the cavities requires an accretion mass, M acc , of 

F 

M acc =^. (3) 

The value of e depends on poorly understood details of the jet 
production process and, probably, on black hole spin. Under 
the usual assumption, that the maximum energy that can be 
extracted is determined by the binding energy of the last stable 
orbit, the upper limit on the efficiency ranges from e < 0.06 
for a nonrotating bl ack hole to e 5, 0-4 for an extreme Kerr 
black hole (Kin2, Frank, & Raine 2002). We assume when 
calculating the energy of the outburst that each cavity repre- 
sents 4pV of energy (i.e. that they are filled with a relativistic 
plasma). 

Since some of the accreting material's mass goes to power 
the jets, the black hole's mass grows by 

AM BH = (l-e)M acc . (4) 
Therefore, increased efficiency results in smaller black hole 
growth for a given outburst energy. The time-averaged ac- 
cretion and black hole growth rates were found by dividing 
Equations | 3| and | 4] by the characteristic time scale discussed 
in Section 13.1.21 Table [3] lists the inferred mass by which 
the black hole grew and the average rate of growth during 
the outburst. The implied black hole growth rates vary across 
our sample by approximately four orders of magnitude, from 
1.6 x 1O" 4 M yr" 1 (M87) to 1.1 M yr" 1 (MS 0735.6+7421), 
with an average value of ~ 0.1 M Q yr" 1 and a median value 
of 0.035 M Q yr" 1 . 

3.3. Eddington and Bondi Accretion Rates 

It is useful to compare the inferred accretion rates to two 
theoretical rates, the Eddington and Bondi accretion rates. 
The Eddington rate is indicative of the maximum likely 
(steady-state) rate of accretion under the assumption of spher- 
ical symmetry, and occurs when the gravitational force acting 
inward on the accreting material is balanced by the outward 
pressure of the radiation emitted by the accretion process. For 
a fully ionized plasma, the Eddington accretion rate is 

M yr-' Z - Ze Wm q J- W 
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TABLE 3 

Black Hole Masses and Growth Rates. 





MBH,meas a 




M BH .L K b 


AM BH C 


Mbh c Bondi ratio 11 Eddington ratio d 


System 


(1O 9 M ) 


(10 9 M S ) 


(10 9 M ) 


(M ) 


(M yr') (M acc /M Bondl ) (Macc/MEdd) 



A85 


, 1+0.6 
1 - 1 -0.4 


1 +0 6 


2 5+2.5 v 


L0 5 


5.9^ ; 9 x Kr 3 


12+52 
lz -9 


2.6!" x 


io- 4 


A133 




7 +0 - 4 
u - '-0.3 


5.0« ; i x 


LO 6 


9.$$% x uy2 


1 000+ 3300 
luuu -670 


7-3+ 2 |g X 


io- 3 


A262 


4 +0 - 2 


0.3 ±0.1 


7 "S +2 -° x 


LO 4 


1.5+i; 2 x Kr 3 


14+41 


7 1+ 3 - 2 x 
Z - 1 -1.0 x 


io- 4 


Perseus 


0.3±0.1 


6 +0 3 
U.D „, 


3.S+ 3 ; 9 x 


LO 6 


2.4+1-5 x 1() -2 


1 400+ 4200 
i+uu_ 900 


3.9+_ 5 ; 2 x 


io- 3 


2A 0335+096 




5+ - 3 


2-2^ x 


L0 5 


3.8+ 3 ] x Kr 3 


36+ 137 
JO -26 


x 


io- 4 


A478 




9 +0 6 


3.0+ 2 ; 3 x 


L0 5 


1.6i>; 2 x Kr 2 


41+206 


9.0+i 7 8 x 


io- 4 


MS 0735.6+7421 




7+ 05 
u -'-0.3 


3.2 +3 -5 x 


10 s 


i i+i.i 
1,1 -0.4 


isooo™ 


7.9+ 2 ° x 


io- 2 


PKS 0745-191 






1.4+" x 


LO 7 


2J+ 2 ; 2 X Kr 1 


630+ii™ 


1-2^ x 


IO" 2 


4C 55.16 




5+ - 3 


2 4+ 2 -5 v 


LO 6 


6.7+™ x ur 2 


1 2 00+ 5600 
lzuu -900 


6.5+i 5 x 


io- 3 


Hydra A 


9 +0 7 


4+ - 2 


1 3+1° v 


LO 7 


6.8+^ x ur 2 


210+j 3 » 


3-711 x 


io- 3 


RBS 797 






7.5t'»x 


10 6 










Zw 2701 




6 +(M 


7.2+|\ x 


10 7 


1 0+ 1 ' 4 
1 - u -0.4 


62000+ 4 «™ 


S-2^0 x 


io- 2 


Zw3146 




2 6+ 2 - 4 


7-7£§ x 


10 7 


9+ 1 1 


370! 3 r 


x 


IO" 2 


A1068 




1 - u -0.4 






3.1 x 10~ 3 


6 4+ 18 


1-6^6 x 


io- 4 


M84 0.36 


0.7 ±0.2 


0.12 ±0.03 


6.1+f l x 


10 2 


1.6+ 2 ; 3 x 10- 4 


44+ 15 
u -++-0.33 


2-11:1 x 


io- 5 


M87 3.3 ±0.7 


1 2+ 05 


0.3 ±0.1 


4 1+ 2 - 9 x 
41 -0.6 X 


10 3 


1.0+* 1 ; 7 x Kr 3 




1.5+»; 7 x 


io- 5 


Centaurus 


0.4 ±0.1 


0.3 ±0.1 


1 7 +1 ° x 


10 4 


1.2+° :9 x ur 3 


2.4!" 


i.e^x 


io- 4 


HCG 62 




0.2 ±0.1 


9 - 3 -5 5 5 X 


10 3 


6.1+ 9 ; 7 x Kr 4 


1 2+ 54 


1.4i 3 ; 7 x 


io- 4 


A1795 


6 +0 3 
u.o_„ 2 




9.4+] 3 X 


10 s 


2.6+ 3 ;* x Kr 2 


oqn+2660 


2.0+f f x 


io- 3 


A1835 






9.5 + i° 3 x 


10 6 


2.8+ 3 ;g X Kr 1 


520^™ 


13 + 4° 7 x 


io- 3 


PKS 1404-267 


0.4±0.1 


0.2±0.1 


2.3 +2 i 9 x 


10 4 


3.2^ x Kr 3 


79+266 
/z -53 


4.3+| : 2 x 


io- 4 


MACS J1423.8+2404 






5.8+J'g x 


10 6 


2.2+|° x Kr 1 








A2029 


, 5+0.8 


2 1+ 1 - 6 


9.6^4 x 


10 s 


1.4+0j x 10- 2 


7 0+ 26 


4.4+«- 3 x 


io- 4 


A2052 


4 +0 - 2 


6 +0 3 

-0.2 


3.6^1 x 


10 s 


2.3+ 3 -J x Kr 2 


510+ 2710 


3.1« : | x 


io- 3 


MKW 3S 




0.3±0.1 




10 6 


6.5+f 7 x Kr 2 


12000+ 98000 
lzuuu -9000 


l.l^x 


io- 2 


A2199 




7 +0 ' 4 
u -'-0.2 


1.5« j x 


10" 


4.3+ 3 ; 9 x Kr 2 


260+ 860 

-170 




io- 3 


Hercules A 




7+ 05 
"•'-0.3 


6.3+ 8 ' x 


10 6 


5.0+f 4 x ur 2 


2100 +i3ioo 

zluu -1600 


3.3!?; 2 x 


io- 3 


3C 388 


1 5+ 12 


6 +0 3 

-0.2 


1-l^x 


10 6 


3.1+| 4 x ur 2 


Q 6 0+ 77 50 


1.0+ 34 x 


io- 3 


3C 401 






2-2^ x 


10 6 


1.0+| 9 x 10"' 








CygnusA 2.7 ±0.7 




1 0+ 06 


1 7+ 1 - 4 x 
1 -'-0.3 X 


10 7 


2.1+ 1 j x 10 -J 


210+630 
zlu -120 


3.6+H x 


io- 3 


Sersic 159/03 




6 +0 4 

-0.2 


5.0+f | x 


10" 


1.2+| 3 x Kr 1 


'+""-1100 


1-0^6 x 


io- 2 


A2597 


0.2±0.1 


0.3±0.1 


7.2+" x 


10 s 


U+ 1 ,- 4 x ur 2 


640+ 4100 


2.6+ 7 ;' x 


io- 3 


A4059 


7+ 10 


1 0+ - 6 


6.1+| 2 x 


10 5 


1.5+i 4 xlO- 2 


450+5830 
^"-410 


1 9+4.6 y 


IO" 3 



a Black hole mass measured using gas kinematics. For Cygnus A, the value of Tadhunter et al. (2003) was adopted, adjusted to our adopted angular diameter- 
distance of 224.2 Mpc. For M87, the average of the values of Harms et al. (1994) and Macchetto et al. (1997) was adopted, adjusted to a distance of 17.9 Mpc. 
For M84, the value of Maciejewski & Binney (2001) was adopted, adjusted to a distance of 15.2 Mpc. 

^Values have been adjusted by a factor of 0.35 (see text for details). 

c The change and rate of change in black hole mass were calculated assuming e =0.1. 

^The Bondi and Eddington rates were calculated with AfBH,meas when available. If no measured value exists, Mbh. ct was used, if available, and Mbh± k if not. 



This rate i s a fu nction only of the black hole mass (discussed 
in Section 13. 4> and the assumed radiative efficiency, e. Ta- 
ble |3] lists the Eddington ratios (M acc /MEdd) for our sample, 
calculated assumin g e = 0. 1 . 

The Bondi rate (Bondi 1952) sets the rate of accretion, as- 
suming spherical symmetry, for a black hole with an accreting 
atmosphere of temperature (T) and density (n e ) as 



MBondi 

M Q yr" 1 



= 0.012 x 



/ kT 
VkeV 



-3/2 



Mi 



BH 



V 10 9 M, 



o 



(6) 



This accretion occurs within the Bondi radius, inside which 
the gas comes under the dominating influence of the black 
hole: 



^Bondi 

kpc 



0.03. x(£ 



Ml 



BH 



V 10 9 M P) 



(7) 



The Bondi rate is therefore an estimate of accretion directly 
from the hot ICM onto the black hole. Table[3] lists the Bondi 
ratios (M acc /MBondi) for our sample, and Table lA2l in the Ap- 
pendix lists the properties used in the calculation of the Bondi 
rates. In calculating the Bondi rate, we use the modeled tem- 
perature and density from Chandra spectra, extracted from a 
central region that contains ~ 3000 counts after the exclusion 
of any non-thermal point sources. However, the size of the 
central region is not sufficiently small to resolve the Bondi 
radius of any system in our sample; therefore, the true tem- 
perature and density of the ICM at the Bondi radius could be 
lower and higher, respectively, than we have measured, re- 
sulting in an underesti mate o f the Bondi rate. We discuss this 
effect further in Sectionl4!2l 



3.4. Black Hole Masses 
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Calculation of both the Eddington and Bondi rates requires 
estimates of the black hole mass. Of the systems in our sam- 
ple, only three (Cygnus A, M84, and M87) have direct mass 
measurements (see Table [3}. For the remaining systems, we 
use the bulge properties of the host galaxy as proxies for the 
black hole mass. As discussed earlier, the black hole's mass 
scales with the large-scale properties of the host galaxy such 
as bulge velocity dispersion and luminosity. The most well- 
studied relation between the black hole mass and the proper- 
ties of the host galaxy is the Mbh - o relation, which relates 
Mbh to the stellar velocity dispersion (a) of the galaxy's bulge 
as 

where a, (3, and Co are constants. The values of these con- 
sta nts vary somewhat fro m study to study (for a discussion, 
see Tremaine et al. 2002). Fo r the purposes of our calcula- 
tions, we adopt the values of [Tremaine et al. (2002), namely 
a = 8.13 ±0.06, 13 = 4.02 ±0.32, and or, = 200 k m s" 1 . 

In deriving this relation, Tremaine et al] J2002I) use as a the 
mean stellar velocity dispersion within a slit aperture of length 
2r e and width 1"— 2" (denoted by <j\). Unfortunately, most of 
our sample lacks dispersions measured in this aperture. In- 
stead, central velocity dispersions (generally measured within 
an aperture of r ~ 2") are more common. Central dispersions 
(denoted by er c ) were taken from the HyperLeda Database. 3 
Measurements of a c exist for 15 of the 33 galaxies in our 
sample (listed in TableQ. When more than one measurement 
exists, we use the weighted average of all available measure- 
ments. We have estimated the magnitude of the error result- 
in g from our use of a c instead of en using the relations given 
in Jdrgensen. Franx. & Kiaergaard ( 1995) and Tremaine et al. 
(2002) and find for the 8 systems in our sample with mea- 
surements of both r e and <r c that Mbh.<t increases on average 
by 10% after the correction, much less than the typical formal 
uncertainties in Mbh, a- Since we lack measurements of r e for 
some systems and the correction is small, we ignore the aper- 
ture correction and use simply a c in our calculation of Mbh.u- 

For the 18 systems without a measurement of velocity 
dispersion, we calculate the black hole mass from the total 
A"-band luminosity of the bulge (Lg) using the relation of 
iMar coni & Hunt (2003) for their group 1 black holes (those 
with secure mass determinations): 



log 



BH,L A - 

M m 



--A + B 



log 



L, 



10.9 



(9) 



where A = 8.21 ± 0.07 and B = 1.13 ± 0.12. Apparent 
^-band magnitudes were taken from the Two Micron All 
Sky Survey (2MASS) catalog. 4 The apparent magnitudes 
were corrected for Galactic e xtinct ion with the values of 
Schlegel. Finkb eTner. & Davisl lll998l) and corrected for red- 
shift (.^-c orrected) and evolution using the corrections of 
PoggiantH dl997l) . Lastly, the magnitudes were converted to 
absolute magnitudes using our assumed cosmology and the 
redshifts listed in Table^] 

We note that there is a systematic offset between the masses 
calculated by the two methods for the 15 systems that have 
measurements of both central velocity dispersion and total K- 
band magnitude. Masses calculated from the total /f-band lu- 
minosity are on average 2.9 ±1.6 times greater. We checked 



3 Available at http://leda.univ-lyon.fr/ 

4 See http://www.ipac.caltech.edu/2mass 



this result using total j?-ba nd m agnitudes from the Hyper- 
Le da database (see Sec tion 13. 6> and the Mbh - Mr relation 
of iMcLure et alJ J200 4) and find a similar sys t ematic offset 
[M BH ,M fi = (3.3 ±2.4) x M BH ,a]. iBettoni et al] J2003) find a 
similar but smaller offset in a sample of radio galaxies and at- 
tribute it to systematically low values of a. Since out values 
of <r are typically weighted averages of several values from a 
number of different sources, it is unlikely that they would be 
systematically low across our entire sample. 

We do not understand the origin of the offset in our data, 
but note that the galaxies in our sample are mostly large cDs, 
with extended stellar envelopes that may bias thei r total mag- 
nitude s with res pect to normal ell i pticals (e.g., IScho mbert 
1986); however, Fuiita & Reiprich (2004) do not find evi- 
dence of such an offset in a similar sample of CDGs. It is also 
possible that the Mrh- c t relation breaks down at high masses 
(see e.g.. IShields et al] 120061) : however, t here is little evi- 
dence to support this hypothesis at this time. Mar coni & Hunt! 
( 2003 ) find evidence of a significant correlation between Mbh 
and the bulge effective radius, with the result that Mbh.o- may 
be too low for large bulges. For typical values of the effective 
radius for galaxies in our sample (r e ~ 10 kpc), the magni- 
tude of th is effect is sufficient to a ccount for the offset we see. 
However, Marconi & Hunt! ( 120 03 ) note that this correlation is 
weak, and further investigation is required to confirm its ex- 
istence. For the purposes of calculating the Eddington and 
Bondi rates, we adjust the black hole masses inferred from 
the K-b&nd luminosities by a factor of 0.35. The black hole 
masses inferred by both methods are listed in Table[3] 

3.5. Star Formation Rates 

The determination of reliable star formation rates requires 
sensitive photometry over a broad wavelength range to iden- 
tify and isolate the star-forming population. Secure star for- 
mation rates are available in the literature for a significant 
number of CDGs. We have collected these rates from the lit- 
erature, adjusted to our assumed cosmology, and their sources 
in Tabled Our sample includes rates derived from both spec- 
troscopic and imaging studies. Readers wis hing to skip the 
technical details should go directly to Section l4~T1 

Typically, in deriving star formation rates, one first finds 
the luminosity of the star-forming population. From broad- 
band images, this luminosity may be found by model- 
ing and subtracting a smooth back ground galaxy (see e.g., 
McNamar a. Wise. & Murravl [20041) . Any extended excess 
emission is then assumed to be due to active star forma- 
tion, and the resulting colors may then be compared to stel- 
lar population models to constrain its age and mass-to-light 
ratio (however, the age and mass-to-light ratio cannot be con- 
strained unambiguously using colors alone). For spectra, a 
similar process is used whereby a spectrum of the background 
galaxy is subtracted (or included as a component in the mod- 
els), and the remaining spectral features are then fit wit h stel- 
lar population models (see e.g.. ICrawford et al] fl999). The 
models constrain the mass-to-light ratio and age of the star 
forming population, which may be used, together with its lu- 
minosity, to calculate the mean star formation rate. In both 
cases, the derived quantities are valid only in the aperture 
used. Consequently, there are three main sources of inho- 
mogeneity in the star formation rates in our sample: differ- 
ences and discrepancies in the model parameters (e.g., as- 
sumed ages), differences in the apertures within which the 
star formation is measured, and uncertainties due to dust ex- 
tinction and reddening. 
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TABLE 4 

Star Formation and Cooling Rates. 



Star Formation Rates Cooling Rates 



Continuous (ref) a Burst (ref) b Aperture XMM RGS (ref) c FUSE (ref) d Aperture 

System (M yr 1 ) (M yr 1 ) (kpc) (M yr 1 ) (M yr 1 ) (kpc) 



A262 


< 0.015 (8) 


r = 0.9 


< 2 (15) 




r = 5 


Perseus 


2.3±0.2 (8) 


r=1.0 




32 ±6 (3) 


11x11 




15.5±5.2 (16) 


r=18.8 










~ 37 (18) 


r = 59 








2A 0335+096 


4.2(17) 


r=16.0 


20±10(15) 




r = 22 


A478 


10.0 (4) 


3.1x44.6 








PKS 0745-191 


16.9±5.6(16) 


r=18.8 








Hydra A 


< 0.5 (7) 


~16(7) r = 4.3 


35 ±20 (15) 




r = 22 


Zw3146 


10.7 (5) 


5.7 x« 25.8) 










< 1 10 (FIR) 










A1068 


18.1 (5) 


3.2x15.6 










28±12 (10) 


46±21 (10) r=10.0 








M84 


< 0.047 (FIR) 






0.32 (2) 


2.2x2.2 


M87 


< 0.02 (8) 


r = 0.2 


< 0.6 (15) 




r = 3.0 




< 0.081 (FIR) 






< 0.44 (2) 


2.6x2.6 


A1795 


0.95±0.10(8) 


r=3.2 




26 ±7 (3) 


36x36 




1.1 (5) 


1.6x7.6 




< 15 (14) 


36 x 36 




2.1±0.9(13) 


19.4x19.4 


< 30(15) 




r=33 




6.3 (9) 


23.2 (9) 13.3 x 26.7 








A1835 


48.9 (5) 


5.1x8.3 


< 200 (15) 




r=99 




79.0(5) 


5.1x8.3 










79.5 (5) 


5.1 x« 23.7) 










140±40(11) 


r = 30 








A2029 


< 0.15 (8) 


r=3.9 




< 27 (3) 


44x44 


A2052 


0.08±0.02 (8) 


r=1.8 


< 10(15) 




r= 17 




0.51 (5) 


0.9x6.8 










0.31±0.10(1) 


r = 2.1 








MKW3S 


< 0.03 (8) 


r = 2.3 


< 10(15) 




r= 11 


A2199 


0.10±0.03 (8) 


r=1.6 










0.10(5) 


0.8x7.3 








A2597 


2.3±1.3 (13) 


16.0x16.0 


<50(12) 




r = 190 




6.4 (9) 


22.3 (9) 24.2 X 24.2 




22 (14) 


48x48 



A4059 



< 10(15) 



r= 147 



REFERENCES. — (1) Blanton et al. 2003; (2) Bregman et al. 2005; (3) Bregman et al. 2006; (4) Cardiel et al. 1998; (5) Crawford et al. 
1999; (6) Lecavelier des Etangs et al. 2004; (7) McNamara 1995; (8) McNamara & O'Connell 1989; (9) McNamara & O'Connell 1993; (10) 
McNamara et al. 2004; (11) McNamara et al. 2006; (12) Morris & Fabian 2005; (13) O'Dea et al. 2004; (14) Oegerle et al. 2001; (15) Peterson 
et al. 2003; (16) Romanishin 1987; (17) Romanishin & Hintzen 1988; (18) Smith et al. 1992. 

a Continuous star formation rate. References are in parentheses. 

^Star formation rate for a burst of star formation calculated as the mass of the burst divided by its age. References are in parentheses. 
c Cooling rates derived from XMM-Newton RGS spectra. References are in parentheses. 
^Cooling rates derived from FUSE spectra. References are in parentheses. 



There are two principle parameters that go into the stel- 
lar population synthesis models: the slope of the initial mass 
function (IMF) and the star formation history. Changes in 
either of these parameters can result in typical deviations of 
factors of ^5-10 in the derived star formation rates. For the 
systems in our sample that have star formation rates avail- 
able, we list in Table 0] rates derived assuming continuous 
star formation for ~ 10 9 yr and, when available, for shorter 
duration bursts. There is little variation across our sample 
in IMF slope, since most studies assume a Salpeter IMF. 
This assumption appea rs to be valid in cooling flows (e.g., 
iMcNamara et al.l2006l) , 

A significant difference between the studies we considered 
is the choice of aperture size. Observations made in spec- 
troscopic slits have the weakness that the star forming re- 
gion may not fall entirely within the slit, resulting in an un- 
derestimate of the total star formation rate. Table @] gives 
the aperture used in each study. In our sample, aperture ef- 
fects could lead to an underestimate of the total star forma- 



tion rate by as much as a factor of ~ 10 if the star formation 
is uniformly distributed acros s the galaxy. However, imag- 
ing studies of CDGs (e.g.. | McNamara & 0'Connei1lll992t 
Cardiel. Gorgas. & Aragon-Salamanca 1998) show that star 
formation is centrally concentrated in most systems, reducing 
somewhat the likely magnitude of this effect. A comparison 
of objects in our sample with star formation rates derived in 
both ways shows that spectroscopic rates are typically lower 
than imaging-derived rates by factors of several. Therefore, 
although spectroscopic estimates should be treated as lower 
limits to the total star formation rates, they are unlikely to be 
more than an order of magnitude lower than the total rates. 

In addition to observational and modeling inhomogeneities 
across our sample, a number of uncertainties exist in any 
derivation of the star formation rate. Principal among 
these are the effects of extinction and reddening due to 
dust. These effects are difficult to quantify without high 
resolution imaging which is not generally available. But 
comparison between the f/-band rates, which are subject to 
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strong extinction, and far IR r ates, which are not, agree to 
within a factor of two (e.g., iMcNamara. Wise. & Murravl 
120041: IMcNamara et all 120061). Of the objects in ou r 
sample. A2052 dBlanton. Sarazin & McNamaral 120031) . 
A1068 (IMcNamara. Wise. & Murravl 12004 . 2A 0335+096 
(Romanishin&Hintzenl 119881) . A1795 and A2597 
(O'D eaetalJ 120041). and those systems studied by 
Crawford et alJ i 19991) have published rates that have 
been corrected for the presence of dust. The intrinsic color 
excess for systems simila r to those in our sample is typically 
E(B-V) ~ 0.3 dCrawford et alll999l) . 

Lastly, errors in mass-to-light ratio and age, while leading 
to errors in accreted mass, generally result in robust star for- 
mation rates due to the compensating effect that older popu- 
lations have higher mass-to-light ratios. Therefore, errors re- 
sulting from an overestimated age will be partly compensated 
by an overestimated population mass, reducing the error in the 
resulting star formation rate. 

A number of objects in our sample have no published opti- 
cal star formation rates, or their rates were measured only in 
small apertures. For these objects, when possible, we have in- 
ferred the star formation rate from the far infrared (FIR) IRAS 
60 /im flux derived with the Infrared Processing and Analy- 
sis Center ' s SCA NPI tool. 5 We used the following relation of 
iKennicutH ( 119981) to convert the total FIR luminosity to a star 
formation rate: 



100.0 



SFR 

M Q yr 



5, 4.5 x 



^FIR 



10 44 



erg s 



(10) 



M\,u\ge = ibulge ( — 



where Lfir ~ LlL^^m (Rowan-Robinson et al. 1997), and 
we have assumed that all the UV photons emitted by young 
stars are absorbed and re-radiated by dust in the FIR. Three 
objects in our sample have reliable 60 /im fluxes: Zw 3146, 
M87, and M84; for these objects we derived upper limits to 
the star formation rates (see Table|4}. 

3.6. Bulge Masses 

Lastly, to estimate the impact of star formation on the mass 
of the galaxy's bulge, we have estimated the mass of the bulge 
as 

) ■ (ID 

bulge 

We use the total /?-band luminosity of the galaxy for the bulge 
luminosity, Lb u i ge - Total apparent magnit udes were taken 
when available from the catalog of Prugn ieT& Heraudeaul 
(1998) and otherwise from the LEDA database (both 
databases are available through HyperLeda). In cases in 
which the /?-band magnitudes were unavailable, we used the 
total B-band magnitudes, if available, and converted these 
to R-b&nd magnitudes using < (B-R) >= 1.44 ±0.17 (the 
average corrected color of our sample). The apparent mag- 
nitudes were corrected and converted to absolute magni- 
tude s in the same way as the A'-band magnitudes (see Sec- 
tion I3.3I >. For the 7?-band mass-to-light ratio of the bulge, 
we adopt {M/L\ n ] m = 6.3 ( M/L)p , the average found by 
iFisher. IlTir igworth. & Franxl ( 119951) . after adjusting to our 
cosmology. The derived absolute magnitudes and bulge 
masses are listed in Tabled 

3.7. Net X-ray and Far UV Cooling Rates 
5 See http://irsa.ipac.caltech.edu/Missions/iras.html 
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FIG. 1 . — Cooling rates derived from XMM-Newton (empty symbols) and 
FUSE (filled symbols) spectra versus those derived from Chandra spectra. 
M87, A1795, and A2597 have b oth XMM-Newton and FUSE rates and hence 
appear twice. M84 is listed in Bregman et al. 1 2005 ) as a probable FUSE 
detection. The dotted line denotes equality between the two rates. 



Gas cooling out of the ICM at T ~ 10 7 K loses its energy 
primarily through thermal emission in the soft X-ray band. 
Therefore, its rate of cooling is best measured from X-ray 
spectra. Grating observations made with XMM-Newton pro- 
vide high spectral resolut ion and hence the best constraints on 
the cooling rates. iPeterson et al.l d2003 ) have derived cooling 
rates from XMM-ftwton grating observations for nine of the 
objects in our sample, and we list the most constraining rate 
(i.e. the smallest rate in any of the temperature bands) from 
this study in Table |4] With the exception of Hydra A and 2 A 
0335+096, these rate s are upper limits . We also list the rate 
derived for A2597 by Morris & Fabian ( 120051) . 

At lower temperatures (T ~ 10 5 K), cooling gas should emit 
strongly in the far ultraviolet, ma inly through the OVI dou- 
blet (see Edgar & Chevalier 1986), where high-quality spec- 
troscopic observations can be made. Such emission has been 
detected by FUSE for an additional six objects in our sam- 
ple. The inferred FUSE cooling rates, calculat ed assumin; 
the OVI emission is due to cooling gas (see e.g., Oegerleetal 
l2001tlBregman et all200H 120061) . are also listed in Table 0] 

We have also derived cooling rate s from lower-spectral- 
resolution Chandra data (see Section l3.1.3> . For compari- 
son, we plot the cooling rates from XMM-Newton and FUSE 
against those from our Chandra analysis in Figure ^ De- 
spite significant differences in aperture and in the details of 
the modeling, the agreement between the X-ray-derived rates 
is reasonably good, as is their agreement with the UV-derived 
FUSE rates. We note, however, that the Chandra rates appear 
to be systematically lower than the other two rates, possibly 
due to spatial and spectral resolution effects or calibration and 
modeling differences. 

It should be emphasized that neither the Chandra nor 
XMM-Newton rates are based on fits to emission lines from 
gas cooling to low temperatures; rather, they are both based 
on fits to the continuum. Additionally, models fit to X-ray 
spectra do not generally require a cooling component to ob- 
tain an adequate fit. Therefore, X-ray-derived cooling rates 
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FIG. 2. — Black hole growth rate versus star formation rate. The diagonal 
lines represent Mbh = 1.4 x W 3 SFR (see text for details) for different values 
of e. Circles denote continuous SFRs measured from broadband images, stars 
denote continuous SFRs measured from spectra taken in slits, and triangles 
denote continuous FIR SFRs. When more than one rate is available, we plot 
the largest rate. If an object has both broadband and spectral rates, we plot 
only the broadband rate. Squares denote rates for a burst of star formation 
and are joined to symbols denoting continuous rates for the same object by 
horizontal lines. 

should be interpreted as the maximum rates of cooling con- 
sistent with the spectra and not as unequivocal detections of 
cooling. Until line emission that is uniquely due to cooling 
below 1 keV is identified, cooling through this temper ature at 
any le vel cannot be confirmed (see however Morris & Fabian 
120051 who find possible weak detections of several cooling 
lines in XMM-Newton data of A2597). However, the reason- 
ably close correspondence between FUSE and X-ray-derived 
rates indicates that cooling is occurring at or just below the 
detection limits. 

4. RESULTS AND DISCUSSION 

4. 1 . Black Holes and Bulges: Simultaneous Growth 

X-ray cavities provide a strong lower limit on the energy 
of the AGN outburst, independent of accretion disk radiation 
models and photon conversion efficiencies. Therefore they 
provide a robust means of estimating the minimum mass ac- 
creted onto the black hole. These properties allow us to inves- 
tigate the relationship between the black hole's growth and 
the bulge's local (small-scale) growth in the same systems in 
a unique and detailed fashion that has not been possible be- 
fore. 

Figure[2]shows the black hole growth rate versus the bulge 
growth rate (traced by star formation) for the systems in 
our sample with reliable star formation rate estimates. We 
plot as dashed lines the time derivative of the present-day 
Magorrian relation, as found by Haring 
1.4 x 10~ 3 Mbuig e . In terms of our derived quantities, this rela- 
tion becomes (1 - e)P cav /(ec 2 ) = 1.4 x IQ~ 3 SFR. 

There is a trend, with large scatter, between the bulge 
and black hole growth rates, centered approximately on the 
Magorrian slope (assuming accretion efficiencies of e ~ 0.1- 
0.4). 

As discussed in Section 13.21 the upper limit on the effi- 



ciency with which the rest mass energy of the accreting ma- 
terial is converted to outburst energy varies between 0.06 
and 0.4. Therefore, under the assumption that star formation 
traces all of the bulge's growth, consistency with general rel- 
ativity and the slope of the Magorrian relation requires that 
all objects with estimates of the total star formation rate fall 
below the e = 0.4 line in Figure|2] 

This requirement is clearly violated in a number of objects. 
For example, the black holes in both Hydra A and A2052 are 
growing faster than strict adherence to the Magorrian relation 
would predict, whereas the bulge of A 1068 is growing too 
fast (note, however, that no cavities were detected in A1068's 
atmosphere; therefore, there is large uncertainty in the rate 
of growth of A1068's black hole). While the discrepancy in 
A1068's rates may be explained with an extremely low effi- 
ciency for the conversion of gravitational binding energy of 
the accreting material to outburst energy (e ~ 0.005), it is also 
possible that present-day growth is occurring in spurts, with 
periods of cooling and star formation (as in A1068) in which 
the bulge grows quickly with little commensurate black hole 
growth, while during periods of heating (as in Hydra A) the 
black hole grows more quickly than the bulge. 

The trend in Figure |2] may be interpreted as an indication 
that, in a time-averaged sense, the growth of the bulges and 
black holes in our sample proceeds roughly along the Magor- 
rian relati on. W hen compared to the bulge masses calculated 
in S ection l3.5l and the black hole masses calculated in Section 
13.21 the black holes are growing at rates of ~ 10~ 9 - 10~ 12 yr" 1 
and the bulges at rates of ~ 10 -11 - 10~ 13 yr" 1 . Present-day 
growth would not move most of the systems significantly off 
of the Magorrian relation, even if growth at such rates was 
constant for the age of the universe. 

However, for a number of systems, current growth could 
produce their present-day black holes in < 10 10 yr. The three 
most extreme cases (MS 0735.6+7421, Zw 2701, and Zw 
3146) have growth rates that, if constant over just ~ 10 9 yr, 
would be sufficient to grow their black holes to their current 
masses. Periodic and powerful outbursts, without commensu- 
rate bulge growth (e.g., MS 0735.6+7421), could cause sig- 
nificant departures from the Magorrian relation. 

These three systems represent ~ 10% of our sample, im- 
plying a duty cycle in active systems of one such outburst 
every ~ 10 8 /0.1 = 10 9 yr. Large outbursts might shut off 
cooling (and hence fuelin g) for long periods, making them 
a relatively rare event (see lDonahue et a l. 2006). If the most 
powerful outbursts are infrequent in the present-day universe, 
the Magorrian relation must have been established during ear- 
lier periods of extreme black hole and bulg e growth, as has 
been postulated by a number of authors (e.g.JYu & Tremairiel 
1 20021 lBinnevll2005tjldTM atteo. Springe l. & Hernquisll 120051 
Churazov et al. 2005) and supported by high-redshift quasar 
studies (e.g. JMcLure & DunlorJ20 04). 

Lastly, it is possible that we are missing some fraction of 
the bulge growth. The CDG is thought to grow through the 
addition of material by tw o main processes: cooling of gas 
out of the ICM (see Section l3.1.1> and merging (cannibalism) 
of the CDG with other cluster member s. The rate of growth 
from mergers is difficult to measure. iLauerl J1988I) . through 
a study of multiple-nucleus CDGs, estimated a cannibalism 
rate of L w 2L* per 5 x 10 9 yr. This estimate implies that such 
growth is significant over the age of the cluster; however, the 
time scale for this growth is much longer than the cooling 
and star formation time scales considered here. Therefore, 
we have neglected mergers and used the star formation rates 
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FIG. 3. — The black hole's relative change in mass versus the mass of the 
black hole, inferred either from gas kinematics or the stellar velocity disper- 
sion (filled symbols) or from the A"-band luminosity of the host galaxy's bulge 
(empty symbols, corrected by a factor of 0.35). 
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FIG. 4. — Bondi ratio (defined as Mace /M-g on <i\) versus the ratio of the semi- 
major axis of the central region (within which the Bondi rate was calculated) 
to the Bondi radius. The symbols are the same as those in Figure |3] The 
lines denote the likely scaling of the measured Bondi ratio with the size of 
the central region, assuming a true Bondi ratio of 1 at the Bondi radius and 
a density profile p oc r , with either a flat core inside o/RBondi = 10 (dashed 
line) or no core (dotted line). 

described in Section l3~5l to set the instantaneous bulge growth 
rate. 

4.2. Accretion Mechanism 

To investigate whether the growth of the black hole, 
or equivalently the energy of its outburst, depends on its 
mass (inferred assuming bulg es come equipped with mature 
black holes, see Section 13. 4> . we plot in Figure [5] the frac- 
tional change in the black hole's mass during the outburst 



(AMbh/A^bh) against its mass. There is no clear indication 
in Figure [5] that the growth of the black hole depends on the 
black hole mass, at least to the extent that the bulge veloc- 
ity dispersion or luminosity is a good black hole mass es- 
timator for these systems. For example, systems that differ 
by a factor of two in inferred black hole mass, such as M84 
(M BH ~ 4 x 10 8 M ) and MS 0735.6+7421 (M BH ~ 7 x 10 8 
Mq), differ in their fractional growth by a factor of ~ 10 5 . 
However, uncertainties in the black hole and accreted masses 
may obscure any underlying correlation. 

For a number of objects in our sample, the implied accre- 
tion rates necessary to generate the cavities are well above our 
Bondi accretion rates (by factors of up to ~ 5 x 10 4 , see Ta- 
bleland Figure^). Specifically, those systems with the most 
powerful outbursts appear to have the largest Bondi ratios, as 
should be expected from the small rang e in black hole masses. 
However, as discussed in Section l3~31 we do not resolve the 
Bondi radius. Therefore, our Bondi rates are probably lower 
that the true values, particularly in higher redshift objects and 
those observations with a low number of total counts (result- 
ing in a larger central region to obtain ~ 3000 counts). 

To illustrate the radial dependence of this correction, we 
plot in Figure |4] the ratio of the accretion to Bondi rate ver- 
sus the semi-major axis of the central region from which the 
Bondi rates were calculated, normalized to the Bondi radius. 
The trend in this figure supports the conclusion that the Bondi 
radius is not resolved. Overplotted are lines denoting the scal- 
ing of the measured Bondi ratio with radius, assuming a true 
Bondi ratio of unity at the Bondi radius and a density profile 
that rises as p oc r~ l to the Bondi radius (upper line) or flat- 
tens inside a/Rncmd\ = 10 (lower line), as observed in M87 
ddi Matteo et alJl2003l) . Objects near or below these lines 
could reasonably have ratios of order unity or less and thus 
be consistent with Bondi accretion. Those significantly above 
the lines are likely to be accreting in excess of their Bondi 
rates. 

All of the objects in our sample are consistent with Bondi 
accretion, but only if the density continues to rise as a pow- 
erlaw to the Bondi radius. The accretion rates in those ob- 
jects with the least powerful outbursts (such as M84 and M87) 
are generally consistent with Bondi ratios of signi ficantly 
less th an unity. This conclusion is supported by Allen et al. 



(2006), who find that accretion rates in ellipticals with low- 
power outbursts are consistent with Bondi accretion. 

However, a number of objects (typically those with pow- 
erful outbursts, such as MS 0735.6+7421 and Zw 2701) are 
barely consistent with Bondi accretion and would have dif- 
ficulty fueling their outbursts through Bondi accretion alone, 
suggesting some other route for much of the accreting ma- 
terial, such as cold accretion (e.g. , the cold feedback mech- 
anism of Pizzolato & Soker 2005). Additionally, the Bondi 
accretion rate assumes spherically symmetric, radial accre- 
tion, while real astrophysical flows will have some net angular 
momentum. An e xample is M87, which appears to posses a 
central disk of gas ( Harm s et alJl994llMacchetto et al.ll997l) : 
thus, any accreting material would be likel y to have significant 
angul ar momentum (for a discussion, see Pizzolato & Sokeil 
2005). Recent hydrodynamic s i mulations of accretion flows 
(e.g., IProga & Begelmanl 120031 iRrumholz. McKee. & Kleii] 
2005 ) find that even small amounts of angular momentum can 
reduce the accretion rate to well below the Bondi rate. 

It is also possible that the Bondi rates (and hence central 
densities) in these objects were higher at the time of the out- 
burst than they are now. We note however that very high den- 
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FIG. 5. — Left: Net cooling rate from XMM-Newton (empty symbols) and FUSE data (filled symbols) versus the star formation rate. Right: Net cooling rate 
from our Chandra X-ray analysis versus the star formation rate. The symbols are the same as those in Figure|2] The diagonal lines denote different ratios of the 
cooling to star formation rate. 



sities imply very short cooling times. At sufficiently high den- 
sities, the gas will cool and fall out of the hot phase in which 
Bondi accretion operates, placing an upper limit on the den- 
sity appropriate f or use in the Bondi calculation (the maximal 
cooling flow, see iNulsen & Fa bian 2000). For example, to 
fuel the outbursts in MS 0735.6+7421 and Zw 2701 by Bondi 
accretion alone, the accretion rate would need to be very close 
to the maximal cooling fl ow value (~ 10% of the Eddington 
rate for these objects, see INulsen & Fabian 2000). However, 
this constraint is not severe enough to rule out Bondi accretion 
as a viable accretion mechanism in most of our sample. 

4.3. Star Formation and Net Cooling of the ICM 

We wish to test the hypothesis that star formation is fueled 
by gas condensing out of the ICM. If true, and cooling and star 
formation vary slowly with time, their rates should be compa- 
rable to each other. To make this comparison, we plot the net 
cooling(condensation) rate against the star formation rate in 
Figure^] Symbols denote the various t ypes of data used to de- 
rive the star formation rate (see Section l33t . The cooling rates 
include estimates inferred from XMM-Newton X-r ay spectra 
and from FUSE ultraviolet spe ctra (se e Section l3~7l . and from 
our Chandra data (see Section l3.1.3> . In almost all cases, the 
X-ray-derived cooling rates should be considered upper lim- 
its, since indisputable evidence of cooling below ~ 1 keV has 
yet to be found in the X-ray emission from cooling flow clus- 
ters. The apparent trend should b e int erpreted cautiously in 
this context (see caveats in Section l3~7l . 

Figure [5] shows the condensation and star formation rates 
for all systems in our sample are consistent with the hypothe- 
sis that star formation is fueled by gas condensing out of the 
ICM. The average ratio of condensation to star formation rate 
for those rates derived in similar apertures is M coo i/SFR ~ 
4, using XMM-Newton and FUSE rates (the ratio does not 
change significantly if Chandra rates are co nsidered) . This 
value is similar to that found by Hicks & Mushotzkv (2005) 
in a study of star formation and cooling using XMM-Newton 
UV monitor data. 



Figure [5] shows that the rates of star formation and cooling 
have converged greatly and are in rough agreement in several 
systems. The classical cooling flow problem, in which the 
X-ray-derived cooling rates were factors of 10-100 in ex- 
cess of the star formation rates in most systems, has largely 
disappeared. While the average discrepancy of four to one 
is still large, it is of the order of the uncertainty in the 
rates. Factors that may contribute to scatter in the rates are 
time-dependent effects such as radio-triggered star formation 
(McNamara & O'Connell 1993) and the time lag required for 
gas at ~ 10 7 K to cool and form stars. 

Lastly, it is clear that if star formation is being fueled by the 
ICM, firm detections of cooling out of the X-ray band should 
be within reach of present and future X-ray obser vatories for 
those objects with large star formation rates ( McNamar a et alJ 
2006). Even with present-day instruments, the upper limits 
on cooling derived to date are very close to the inferred total 
star formation rates for a number of objects (e.g., A1835). 
If this star formation scenario is to survive, future deep X- 
ray observations of these objects (with Constellation-X, for 
example) should detect this cooling gas. 

4.4. Quenching Cooling Flows 

We have demonstrated that, in many systems, the net ICM 
cooling (condensation) rate is in rough agreement with total 
star formation rate. However, we have not dealt with the ques- 
tion of what maintains the bulk of the ICM at X-ray temper- 
atures, preventing it from cooling out at the expected classi- 
cal rates (typically ~ 10- 100 times the star formation rates). 
AGN heating, through cavities, shocks, and sound waves, has 
emerged as the favored mechanism to prevent this massive 
cooling in cooling flows. To investigate whether AGN cavi- 
ties are powerful enough to balance the radiation emitted by 
the ICM, we plot in Figure [6] the cavity power of the central 
AGN against the total radiative luminosity of the intracluster 
gas within the cooling radius (minus the luminosity due to net 
cooling, given in T able Et- This plot supersedes Figure 2 of 
iBirzan et alJ J2004I) . to whose sample we have added deeper 
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FIG. 6. — Cavity power of the central AGN versus the X-ray luminosity 
of the intracluster medium inside the cooling radius that must be offset to 
be consistent with the spectra (Licm = Lx~ L coo \). The symbols and wide 
error bars denote the values of cavity power calculated using the buoyancy 
timescale. The short and medium-width error bars denote the upper and 
lower limits of the cavity power calculated using the sound speed and refill 
timescales, respectively. Different symbols denote different figures of merit: 
circle - well-defined cavity with bright rims, triangle - well-defined cavity 
without bright rims, square - poorly defined cavity. The diagonal lines denote 
Pcav = ^icm assuming pV, 4pV, or Y6pV as the total enthalpy of the cavities. 



X-ray data and 14 new cavity systems, most of which lie in 
the upper half of cavity powers. 

Remarkably, most of the systems in our sample have cav- 
ity powers sufficient or nearly sufficient to balance the en- 
tire radiative losses of the ICM within the cooling radius. 
The remaining systems may require other forms of heat 
to offset cooling com pletely, such as thermal conduction 
JVoigt & Fabianll200l . However, we note that the time- 
dependent nature of AGN feedback does not require that cool- 
ing is always balanced by heating. It is possible that those 
systems that do not currently balance are in a cooling phase 
and will be entering a heating phase soon. Intermittent heat- 
ing and cooling would allow for cooling and star formation at 
observed levels. 

In Figure |6] a number of systems lie well above the ApV 
line and even above the pV line, implying that their cavities 
likely represent more energy than required to balance cool- 
ing. These systems, many of which possess supercavities and 
shocks extending beyond the cooling radius, have enough en- 
ergy to quench cooling and to contribute to cluster preheat- 
ing. An example is MS 0735.6+7421, the most powerful such 
outburst know to date. The AGN in this cluster has dumped 
~ 1/3 keV per particle into the IC M (including the energy of 
the shock; McNamar a et al.l l2005). The cavities alone have 
enough energy to quench cooling 15 times over. This amount 
of energy, even if distributed only partly inside the cooling ra- 
dius, should have a profound effect on any cooling gas. Such 
objects may thus be in a heating phase. 

However, unresolved problems still remain, the most obvi- 
ous of which is the absence of cavities and star formation in 
many cooling flow systems. The absence of cavities currently 
does not however rule out significant feedback in the past. 



iDonahue et alJ J2006I) find elevated entropy levels throughout 
the cooling region of both cooling-flow clusters that show evi- 
dence of AGN feedback and those that do not, consistent with 
a history of AGN feedback in all cooling flows. Secondly, 
in cooling flows that lack evidence of AGN heating, other 
sources of heat, such as thermal conduction, may be impor- 
tant. Lastly, the most powerful explosions, such as seen in 
MS 0735.6+7421, may turn off accretion and hence AGN ac- 
tivity for extended periods. 

5. CONCLUSIONS 

We have presented an analysis of the star formation and 
AGN properties in 33 CDGs in the cores of cooling flows. 
We find that the AGN outbursts in most of the systems have 
enough energy to offset most of the radiative losses of the 
ICM, and to severely reduce cooling to levels that are ap- 
proaching the star formation rates in the central galaxy. Using 
the cavities to infer black hole growth and star formation to 
infer bulge growth, we find that bulge and black hole growth 
rates scale roughly with each other in rough accordance with 
the slope of the Magorrian relation. The large scatter may in- 
dicate that growth occurs in spurts, with periods of cooling 
and star formation interspaced with periods of heating, or that 
the efficiency of the conversion of the gravitational binding 
energy of the accreting matter to outburst energy varies across 
the sample. We find the central supermassive black holes are 
growing at rates of ~ 10~ 4 M yr" 1 to ~ 1 Mq yr" 1 (with a 
median rate of 0.035 M© yr -1 ), which, in most of our sample, 
are insufficient to account for their present-day masses. How- 
ever, a number of black holes are growing at rates that are 
consistent with their formation from scratch in ~ 10 10 yr. The 
extreme cases are those objects experiencing the most power- 
ful outbursts (P cav ~ 5 x 10 45 ergs s" 1 , approximately 10% of 
our sample), which are growing at rates sufficient to assemble 
their black holes in ~ 10 9 yr. 

Across our sample, the inferred black hole accretion rates 
are well below their Eddington limits but above their Bondi 
rates. Chandra does not resolve the Bondi radius in these 
systems, and thus significant Bondi accretion cannot be ruled 
out. The exceptions are those systems with powerful out- 
bursts, where either direct accretion from the hot ICM is not 
the principle route of cooling gas or their central properties 
were very different at the time of the outburst than those of 
typical nearby CDGs such as M87. 

We test the scenario that the active star formation is fueled 
by cooling (condensation) from the ICM. We find that star for- 
mation and cooling rates are converging (to an average ratio 
of cooling to total star formation rate of four to one), and in 
some cases are consistent with one another. Inhomogeneities 
in star formation rates and the lack of firm detections of cool- 
ing in X-ray data are the main factors that limit our conclu- 
sions. Nevertheless, this rough agreement is far different from 
the situation a decade ago, when the best X-ray cooling rates 
were tens to hundreds of times the star formation rates. 

Using the best X-ray data to da te, we extend and revise 
the heating versus cooling plot of Birzan et al. (2004). Re- 
markably, we find that AGN heating, as traced by the power 
in X-ray cavities alone, is capable of balancing the radiative 
losses of the ICM in more than half of the systems in our 
sample. However, the means by which the AGN's jet energy 
is converted to heat in the ICM and the efficiency of this con- 
versio n are not yet clear (e.g.. lRevnolds. Heinz. & Begelman 
2002). Additionally, our estimate of AGN heating neglects 
other significant sources of heat that are likely to be present in 
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many of the systems in our sample, such as weak shocks (e.g., 
McNamai-aiSD IZQQj iNulsen et al.ll2Q05albt iForman et alJ 
2005), sound wa ves (e.g.. | Fab ian et al.1 l2006l) . and thermal 
conduction (e.g., Voigt & Fabian 2004). All of these heat 
sources may play a role in maintaining the rough balance of 
heating to cooling observed to exist throughout the cooling 
region. AGN, however, have emerged as the most important 
heating mechanism in cooling flows. 

A unified picture of star formation, cooling, and AGN feed- 
back is now emerging, one with applications to the more gen- 
eral problems of galaxy formation and the truncation of the 
high end of the luminosity function of galaxie s. Both sim- 
ulatio ns and models of galaxy f o rmation (e.g..lBalogh etalJ 
1200 11 ISiiacki & Springell 120061 IVoit. Kav. & Brvanl I2005F 
conclude that AGN heating is required to prevent the over- 
cooling problem in CDM models, in which too many large 
galaxies are formed. A plausible scenario is that AGN regu- 
late the cooling of gas in the cores of cooling flows, preventing 
most of the ICM from cooling but allowing some net conden- 



sation that feeds both star formation and black hole growth, 
possibly in an intermittent manner, along the Magorrian rela- 
tion. 

In summary, we find that star formation and cooling rates 
and AGN outburst energies for our sample of CDGs are 
broadly consistent with the simple AGN-ICM feedback sce- 
nario in which gas cooling out of the ICM feeds AGN out- 
bursts that heat the gas in the cluster's core. Some low-level, 
net cooling may still proceed, and upper limits on its rate are 
consistent with the scenario in which net cooling is the source 
of material for active star formation, unusual in most ellipti- 
cals but present in many CDGs. 
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APPENDIX 

CAVITY AND CENTRAL ICM PROPERTIES. 

Table lAlllists the properties of each cavity measured from Chandra images. Errors in pV include an estimate of the projection 
effects (see Birz an et al.l2004l) . 



TABLE A 1 
Cavity Properties. 





Cavity 


a b 


b c 




pV 


K 


fretill 


^buoy 


System 


FOM a 


(kpc) 


(kpc) 


(kpc) 


(10 58 erg) 


(10 7 yr) 


(10 7 yr) 


(10 7 yr) 


A85 


2 


8.9 


6.3 


21 


1 2 +L2 


2.3 


5.1 


4.2 


A133 


3 


41 


21 


32 


24+11 


3.8 


14 


5.1 


A262 


2 


5.4 


3.4 


8.7 


n 060+° 050 
UUDU -0.017 


1.5 


2.9 


1.7 




2 


5.7 


3.4 


8.1 


n 065+ 048 


1.4 


2.8 


1.6 


Perseus 


1 


9.1 


7.3 


9.4 


3.7+tl 


1.0 


4.9 


1.6 




1 


8.2 


4.7 


6.5 


1 6+ L0 


0.7 


3.6 


1.1 




2 


17 


7.3 


28 


3 9+ 3 - 5 


3.1 


10 


8.3 




2 


17 


13 


39 


9.7+'°/ 


4.0 


13 


10 


2A 0335+096 


2 


9.3 


6.5 


23 


1 0+''° 

1 - u -0.3 


3.2 


6.3 


5.4 




3 


4.8 


2.6 


28 


072+ - 037 


3.7 


4.6 


11 


A478 


2 


5.5 


3.4 


9.0 


74+ - 57 


1.0 


2.9 


1.8 




2 


5.6 


3.4 


9.0 


76+° 55 


1.0 


3.0 


1.8 


MS 0735.6+7421 


2 


110 


87 


160 


77Pt+960 
' /u -360 


13 


58 


26 




2 


130 


89 


180 


830+™ 


15 


66 


33 


PKS 0745-191 


3 


26 


17 


31 


69+* 


3.0 


12 


5.2 


4C 55.16 


2 


10 


7.5 


16 


4.7+« 


1.7 


5.6 


3.0 




2 


13 


9.4 


22 


7.1+™ 


2.3 


7.3 


4.1 


Hydra A c 


2 


18 


12 


29 




3.0 


8.7 


5.1 




2 


20 


12 


31 




3.2 


9.3 


5.6 




3 


42 


21 


78 


20+ 8 


7.8 


21 


17 




3 


34 


24 


66 


27+27 


6.6 


19 


12 


RBS 797 


2 


13 


8.5 


24 


18+' 4 


2.2 


7.5 


5.0 




2 


9.7 


9.7 


20 


2 ( )+36 
zu -13 


1.8 


6.5 


3.4 


Zw 2701 


2 


46 


41 


54 


220+ 340 
zzu -130 


5.2 


23 


7.8 




2 


39 


34 


49 


1 30+ 190 


4.7 


20 


7.2 


Zw3146 


2 


51 


21 


40 


1 70 +i so 


3.7 


17 


6.8 




2 


36 


30 


59 


210 +280 


5.0 


21 


10 


M84 


2 


1.6 


1.6 


2.3 


002+ 004 
u - uuz -0.0015 


0.5 


1.0 


0.4 




2 


2.1 


1.2 


2.5 


oo 1+0.001 
u - uul -0.0005 


0.6 


1.0 


0.5 


M87 


2 


2.3 


1.4 


2.8 


OOifi+0012 
uulo -0.003 


0.4 


0.9 


0.4 
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System 


Cavity 
FOM a 


(kpc) 


b c 
(kpc) 


(kpc) 


pV 
(10 58 erg) 


(10 7 yr) 


frcflll 

(10 7 yr) 


'buoy 

(10 7 yr) 




2 


1.6 


0.8 


2.2 


0.004-0-002 





.4 


0.7 


0.4 


Centaurus 


1 


3.3 


2.4 


6.0 


0.038+0™ 


1 


.0 


2.2 


1.3 




1 


3.3 


1.6 


3.5 


022+ - 012 
u - uzz -0.003 





.6 


1.5 


0.7 


HCG 62 


2 


5.0 


4.3 


8.4 


u - uz '-0.015 


1 


.8 


2.9 


1.5 




2 


4.0 


4.0 


8.6 


() rnn+0.034 
uuly -0.013 


1 


.8 


2.8 


1.6 


A1795 


3 


19 


7.2 


19 


4.r« 


1 


.9 


6.8 


3.7 


A1835 


3 


16 


12 


23 


27.58 


2 


.1 


8.3 


4.1 




3 


14 


9.7 


17 


20+20 


1 


.5 


6.5 


2.7 


PKS 1404-267 


2 


3.5 


2.6 


4.6 


054+ - 060 





.8 


2.0 


0.9 




2 


3.2 


2.7 


3.8 


062+ 0085 
uuoz -0.031 





.6 


1.8 


0.6 


MACS J1423.8+2404 


2 


9.4 


9.4 


16 


15 +27 


1 


.5 


5.7 


2.5 




2 


9.4 


9.4 


17 




1 


.6 


5.9 


2.8 


A2029 


3 


13 


7.2 


32 




2 


.5 


6.8 


6.9 


A2052 


1 


11 


7.9 


11 


1 2+ L4 


1 


.8 


5.5 


1.9 




1 


6.5 


6.2 


6.7 


53+°- 88 


1 


.2 


3.6 


1.0 


MKW3S 


3 


54 


23 


59 


38+f 


6 


.0 


21 


12 


A2199 


2 


15 


10 


19 


3.7+3-7 


2 


.1 


7.1 


3.2 




2 


16 


10 


21 




2 


.3 


7.7 


3.8 


Hercules A 


3 


26 


21 


60 


13+' 8 


6 


.1 


18 


13 




3 


47 


19 


58 


18+| 2 


6 


.0 


19 


13 


3C388 


2 


15 


15 


27 


2 9+ 53 
z - y -1.9 


2 


.9 


7.6 


3.6 




2 


24 


10 


21 


2 3+ 2 - 2 


2 


.4 


6.9 


3.1 


3C401 


2 


12 


12 


15 


5.4+9J 


1 


.6 


6.4 


2.1 




2 


12 


12 


15 


5.4+| 8 


1 


.6 


6.4 


2.1 


Cygnus A 


1 


29 


17 


43 


28+J 8 


3 


.4 


15 


8.5 




1 


34 


23 


45 


56+« 


3 


.6 


17 


7.8 


Sersic 159/03 


3 


20 


14 


23 


10+| 


2 


.9 


9.3 


3.8 




3 


22 


17 


26 


15+ 17 


3 


.3 


11 


4.2 


A2597 


2 


7.1 


7.1 


23 


1 5+ 26 
1,J -0.9 


2 


.5 


7.9 


6.8 




2 


10 


7.1 


23 


2 1+ 20 


2 


.4 


8.6 


6.6 


A4059 


2 


20 


10 


23 


9 o+l.O 
z - z -0.3 


2 


.7 


8.4 


4.2 




2 


9.2 


9.2 


19 


0.84+ 1 ;* 


2 


.3 


6.2 


3.5 



a Figure of merit. The FOM gives a relative measure of the cavity's contrast to its surroundings: (1) high contrast: bright rim surrounds cavity; (2) medium contrast: bright rim 
partially surrounds cavity; and (3) low contrast: no rim, or faint rim surrounds cavity. 
Projected semi-major axis of the cavity. 
c Projected semi-minor axis of the cavity. 

Projected radial distance from the core to the cavity's center. 
e The deeper image of Wise et al. (2006, in preparation) of Hydra A shows two large outer cavities beyond those measured here, but for consistency we report only those cavities 
apparent in archival data. 

Table IA2I lists the modeled temperature and density of the central region (with semi-major axis, a) used in the calcula- 
tion of the Bondi rate. Also listed is the Bondi radius, calculated (using Equation^ from the central temperature and the black 
hole mass given in Table[3] 

TABLE A2 



Central ICM Properties. 




kT 


n e 


a 




System 


(keV) 


(cm" 3 ) 


(kpc) 


(kpc) 


A85 

A133 

A262 

Perseus 

2A 0335+096 

A478 


0.86+00,1 

4 4+0.5 
^-0.4 
I 4+0.1 

9 7+0.3 
z - '-0.3 


o.io7!Sg 

0.048+0-004 
065+° 008 

uut "-0.007 

ulJU -0.005 
n 056+0 .003 

90+0.01 
u - zu -0.02 


5.8 
8.0 
3.4 
8.6 
5.1 
5.3 


0.017 
0.012 
0.013 
0.004 
0.012 
0.010 
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TABLE A2 — Continued 





kT 


«e 


a 




System 


(keV) 


(cm 4 ) 


(kpc) 


(kpc) 


"R T C 1 All C /~ i —7 A O 1 

MS 0735.6+7421 


J -0.2 


n nt-7+0 002 

0.067;"™ 2 


23.8 


0.007 


FKS 0/45-191 


o £+0 4 
2 - 6 -0.4 


A 1/1+001 

0.14_ om 


11.2 


0.013 


Hydra A 


2 - 6 as 


A 1 C+0.0I 

U ' i:, -0.02 


4.7 


0.01 1 


Zw 2701 


T 1+0 3 

* S " 3 -0.3 


n no ,1+0 002 

°- 024 :o:oo2 


37.6 


0.006 


Zw 3146 


t 1+0.3 
J ' I -0.2 


n i -7-7+0 007 
°" 177 -0.007 


15.0 


0.027 


M84 


n en+0 01 

°- 57 ±ao! 


r\ i r»c+0 007 

°- 105 -o:oo7 


0.9 


0.020 


M87 


n n/i+0 02 


r> mi+0 009 
u,±71 -0.009 


1.0 


0.110 


Centaurus 


n -7-7+0 01 

°- 77 -aoi 


A T3+0 0I 


1.3 


0.015 


HCG 62 


n £-7+0 01 


A AC-7+0 007 

0.057;";^ 


2. 1 


0.010 


A 1795 


n -7+0 6 
2J OA 


A At-7+0 00 s ) 

°- 067 -o:oo5 


9.5 


0.007 


A1835 


A A+0 3 

4 -° 03 


r> i 1 /-i+O 003 

°- 11( to:oo3 


27.2 


0.015 


PKS 1404-267 




°- 046 -aool 


8.5 


0.009 


A2029 


o n+0 3 

2 - 9 ±a2 


A 1-7+0 04 

037 -o:o3 


2.2 


0.022 


A2052 


n -71+0.04 


n a 1 -7+0 002 
0.017_ 0002 


5.5 


0.017 


MKW 3S 


9 q+0.8 
2 - S -0.5 


n mo+0 006 
28 -0.009 


7.8 


0.003 


A2199 


9 9+0.2 
z ' z -0.l 


099+° 005 
ULW -0.005 


4.4 


0.010 


Hercules A 


2 +0 ' 2 


0111 +0.0006 
u.uill 00005 


67.0 


0.012 


3C388 


3 O 4 *- 2 

■'•"-O 2 


0069+° 0004 
u ' u -0.0004 


55.6 


0.016 


Cygnus A 


5 2^-5 
-'• z -0.6 


n n9+0.009 
u - 1JZ -O.008 


5.3 


0.017 


Sersic 159/03 


°-0.I 




12.2 


0.010 


A2597 


1,D -0.2 


n 07^+0.005 
U - U/J -O.005 


11.0 


0.006 


A4059 


2 l +ai 


022+ 001 


10.6 


0.010 
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